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A proposal to detect the purported canted antiferromagnet order for the v = 0 quantum Hall state of graphene 
based on a two-terminal spin transport setup is theoretically discussed. In the presence of a magnetic field 
normal to the graphene plane, a dynamic and inhomogeneous texture of the Neel vector lying within the plane 
should mediate (nearly dissipationless) superfluid transport of spin angular momentum polarized along the z 
axis, which could serve as a strong support for the canted antiferromagnet scenario. Spin injection and detection 
can be achieved by coupling two spin-polarized edge channels of the |v| = 2 quantum Hall state on two opposite 
ends of the v = 0 region. A simple kinetic theory and Onsager reciprocity are invoked to model the spin 
injection and detection processes, and the transport of spin through the antiferromagnet is accounted for using 
the Landau-Lifshitz-Gilbert phenomenology. 


Introduction .—Unique electronic properties of 
graphene (a monolayer of graphitic carbon) stem from its 
hexagonal lattice structure, engendering relativistic effects 
at electronic velocities well below the speed of light (T| . 
Graphene is the thinnest and the strongest of 2D materials, 
and an outstanding electrical and heat conductor, holding 
great promise as a building block for future electronic 
devices m. 

A hallmark of graphene’s electronic properties is man¬ 
ifested in magnetotransport. The integer quantum Hall 
(QH) sequence with anomalous filling fractions v = ±4 (n+ 
i/2) mm directly reflects the weakly-interacting massless 
relativistic (“Dirac”) nature of its low-energy excitations 
and the fourfold degeneracy associated with the electron 
spin and valley isospin. The valley degree of freedom dis¬ 
tinguishes between the two inequivalent “Dirac points” in 
the Brillouin zone where the conduction and valence bands 
of graphene touch Bl . 

Under high magnetic fields, strongly-correlated QH 
phases can emerge EIE including the v = 0 state at the 
charge neutrality point, which indicates that interaction- 
induced SU(4)-symmetry breaking within the spin-valley 
space lifts the fourfold degeneracy of the zeroth Landau 
level iHia. A challenge is to understand precisely how this 
symmetry is broken for the v = 0 state. Charge-transport 
experiments, utilizing both the two-terminal and Hall-bar 
geometries, suggest that the bulk and edge charge excita¬ 
tions for the state are gapped 0. Furthermore, the recent 
observation of gapless edge-state reconstruction in tilted 
magnetic field 0 is consistent with the scenario where the 
v = 0 ground state is a canted antiferromagnetic (CAF) 
insulator ED. where the spins Sa on sublattice A have a 
different orientation relative to the spins S B on sublattice 
B\ in the presence of an external magnetic field B normal 
to the graphene plane (defined to be the xy plane), the total 
spin Sa + S B points antiparallel to the field while the Neel 
vector Sa~S b lies in the graphene plane. Despite these re¬ 
cent developments, a more direct experimental verification 


of this CAF scenario would be highly desirable. 

Essentially disjoint from the field of graphene QH 
physics, spintronics is witnessing an increasing interest in 
realizing spin transport through magnetic insulators via co¬ 
herent collective magnetic excitations, which allow for su¬ 
perfluid (nearly dissipationless) transport of spin OllO). 
A recent theoretical work has shown that such superfluid 
spin transport can be realized in antiferromagnetic insula¬ 
tors using a two-terminal setup Do): by laterally attach¬ 
ing two strongly spin-orbit-coupled normal metals at two 
opposite ends of the insulator, both spin injection and de¬ 
tection could be achieved via electrical means using direct 
and inverse spin Hall effects. Transplanting this idea to the 
purported v = 0 CAF state in graphene, superfluid trans¬ 
port of spin polarized along the z axis could be harnessed 
by the CAF via a dynamic Neel texture that rotates about 
the z axis within the graphene plane not . The observa¬ 
tion of such spin superfluidity in the v = 0 state should, 
therefore, support the above CAF scenario. 

Superfluid spin transport .—To realize and detect super¬ 
fluid spin transport through the v = 0 CAF we propose 
the device shown in Fig. |TJa), where the central CAF re¬ 
gion is sandwiched by two v = -2 QH regions; we ig¬ 
nore the effects of thermal fluctuations of the spins in the 
CAF. Spin injection into the CAF is achieved using the two 
co-propagating edge channels of the left v = -2 region. 
Based on the theory of QH ferromagnetism ifTTl these 
edge channels, away from the injection region (shaded in 
green), are in oppositely polarized spin states (labeled f,!) 
collinear with the external field (along the z axis) [the in¬ 
jection region contains the vertices , two corners of the edge 
channels labeled a and b (represented by red circles in 
Fig. [lja)), and the line junction , the stretch of edge chan¬ 
nels between the vertices, adjacent to the CAF]. In the 
absence of spin-flip processes, the edge channels undergo 
very little equilibration outside the injection region fT2l SO 
that their voltages, as they enter the region, are defined by 
the reservoirs from which they originate, i.e., V^. A possi- 
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FIG. 1. (a) Theoretically proposed setup for realizing and de¬ 
tecting spin superfluid transport through the v = 0 QH state of 
graphene. Boundaries between regions with different filling frac¬ 
tions are represented by dashed lines. Two independently biased 
spin-polarized channels on two opposite sides of the v = 0 re¬ 
gion are used to inject and detect spin current flowing through 
the y = 0 region. The external field B constrains the Neel vec¬ 
tor to lie within the xy plane. A slight misalignment of the edge 
spins away from the z axis in the injection and detection regions 
may result from the neighboring CAF (see the blow-up illustra¬ 
tion). The spin states of the v = -2 edge channels are polarized 
collinearly to the z axis outside the injection and detection re¬ 
gions. (b) One possible experimental realization of the proposed 
setup in (a). Both spin channels biased at impinge upon the 
gated v = -1 region, which filters one of the spin channels (the 
spin-down channel) from entering the inner v = -2 region. The 
spin-down channel within the inner region can be separately bi¬ 
ased by Vi thereby allowing independent control of the electro¬ 
chemical potentials of the two spin channels. An analogous setup 
is used for the detection side. 

ble experimental realization of the proposed setup is shown 
in Fig. |TJb), which also shows how the two spin channels 
can be separately biased. Both spin channels biased at Vf 
impinge upon the gated v = -1 region, which filters one 
of the spin channels (the spin-down channel) from entering 
the inner v = -2 region. The spin-down channel within the 
inner region can be separately biased by V|, allowing in¬ 
dependent control of the electrochemical potentials of the 
two spin channels that impinge upon the injection region. 
An analogous setup can be used for the detection side. 

When Vf > V|, inter-channel scattering may occur in- 
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FIG. 2. A cartoon of the CAF in the dc superfluid state. The 
Neel vector rotates within the graphene plane about the z axis 
with a global precession frequency Q, and harnesses superfluid 
spin transport. The phase (the azimuth) of the Neel vector may be 
spatially inhomogeneous if spin current is non-uniformly injected 
along the injection region. Spin current injected via the injection 
region has a static contribution I s oc V_ = V T - and a dynamic 
(spin-pumping) contribution i s oc fl that pumps spin current back 
out into the edge. Two analogous contributions exist also on the 
detection side. 


side the injection region, entailing redistributed charge cur¬ 
rents, 7y and I\, emanating from the region and a net loss of 
spin angular momentum, polarized along the z axis, inside 
the region. Neglecting any external sources of spin loss 
(e.g., spin-orbit coupling, magnetic impurities, etc.) inside 
the injection region, the net spin lost in the edge should 
be fully absorbed by the CAF, leading to the injection of 
spin current (hereafter always defined to be the component 
polarized along the z axis) into the CAF. This will eventu¬ 
ally induce the CAF into a dynamic steady-state, in which 
the local Neel vector in the CAF precesses about the z axis 
with a global frequency O (see Fig. 13 Go). The dynamic 
Neel texture, in turn, pumps spin current lfl3l out into the 
detection edge channels in the detection region, thus facil¬ 
itating the superfluid spin transport from the injection to 
the detection side (the detection region, involving vertices 
a' and b', is shaded in blue in Fig. [lja) and defined anal¬ 
ogously to the injection region). Away from the detection 
region, the spins of the edge channels in the right v = -2 
region are also oppositely polarized along the z axis, so 
that the spin current ejected into the detection edge can be 
determined by measuring the difference in spin current en¬ 
tering and exiting the detection region. 

Theory and results. —A theory for superfluid spin trans¬ 
port through antiferromagnetic insulators has been devel¬ 
oped in detail in an earlier work by the authors Go), and 
the discussion below follows directly from the work. Once 
the dynamic steady-state is established in the CAF, the to¬ 
tal spin current I s entering the CAF via the injection region 
has two contributions: I s - I s + i s , where I s is the spin 
current injected into a static CAF in equilibrium, and i s is 
the spin-pumping (dynamic) contribution describing spin 
current pumped back out to the edge due to the nonequi¬ 
librium Neel dynamics (see Fig. [2]) ff3l . Within linear- 
response, once the static contribution is quantified, the dy- 
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namic contribution can be determined using Onsager reci¬ 
procity, as we show below. 

The static contribution to the spin current within linear- 
response reads I s = (h/2e)[g Q (V-f - V x ) - (/ T - /[)], 
where gQ = e 2 /h and e > 0 is the magnitude of the 
electron charge; 7^ denotes the charge currents emanat¬ 
ing from the injection region in the static limit. Due 
to charge conservation, and the fact that equally-biased 
edge channels leads to equal outgoing charge currents (i.e., 
V T = V| implying f = 7^), the charge currents emanat¬ 
ing from the injection region can be written generally as 
la = gQ[V + + <x(l - y)V.\/ 2, where V ± = V T ± V; and 
cr = ± corresponds to the T and | channels, respectively. 
The real parameter 0 < y < 1 characterizes the strength of 
inter-channel scattering in the injection region (it is explic¬ 
itly computed using a simple microscopic model at a later 
point in this work). The limit of no inter-channel scatter¬ 
ing corresponds to y = 0, while the limit of strong scatter¬ 
ing (full equilibration between the channels) corresponds 
to y - 1. Inserting l a into the expression for I s , one ob¬ 
tains 


l = Ye gQyV - ■ (1) 

The dynamical contribution i s follows from Eq. Q and 
Onsager reciprocity. Let us first define two continuum 
variables in the CAF that are slowly varying on the scale 
of the magnetic length: n{x) and m(x ), n(x) being a unit 
vector pointing along the local Neel order and m(x) be¬ 
ing the local spin density. The global frequency O of the 
rotating Neel texture effectively acts as an additional mag¬ 
netic field in the z direction and introduces a uniform fer¬ 
romagnetic canting of the CAF spins along the z direction 
in addition to the existing canting due to the external field. 
Therefore, in the dynamic steady-state the CAF is char¬ 
acterized by a uniform m{x) = m z e z . Defining the total 
spin M z = m z LW , where L and W are the dimensions of 
the CAF region [see Fig.[lja)], the dynamics of M z in the 
presence of the injected spin current I s is given by 

M Z = I S + •••, ( 2 ) 

where the ellipsis denotes terms arising from the intrinsic 
dynamics within the CAF. Inserting the static contribution 
Eq.Q in for I s in Eq. ([2]) introduces terms linear in V j and 
V|, which are the forces conjugate to the charge currents f 
and 7|, respectively. Onsager reciprocity then endows the 
static contributions l a with a dynamic contribution as 

fi 

I<r = la- V — gQyfviz > ( 3 ) 


It) 



FIG. 3. A cartoon energy diagram at a v = 0 to v = -2 transi¬ 
tion region. In the v = -2 region, energies of the two spin states, 
oppositely polarized along the z axis, are drawn; the Zeeman ef¬ 
fect gives an energetic advantage to the spin-down state. In the 
v — 0 region, the two occupied branches of the CAF spectrum 
are shown. There, an external field in the positive z direction 
results in a ferromagnetic canting of spins in the negative z di¬ 
rection inside the antiferromagnet. Spin orientations of the chiral 
modes are intermediate between the up and down spin eigenstates 
within the v = -2 region (left side) and the canted spins within 
the CAF (right side). The black lines are merely a rough guide for 
the energies of the spin states in the transition region. The above 
illustration does not contain two other branches of the spectrum 
that are a part of the zLL but not essential for the edge physics in 
the transition region. 


axis]. Noting that the force /m z relates to the local Neel 
vector via f Mz = -(n x fi) • e z « -Q da, the total injected 
spin current I s = (h/2e)[gQ(V^ - Vf) - - /|)] can be 

obtained using Eq. ^ as 

I s = f(eV- - Ml) = I, + i,. (4) 

Based on a fully analogous consideration on the detection 
side, the total spin current injected into the edge from the 
CAF becomes /' = ~(y' /4n)(eV'_ - fiQ) = 7' + , where y' 
is the inter-channel scattering parameter, analogous to y, 
for the detection side. Fixing the voltages of the electron 
reservoirs on the detection side to zero, i.e., V'^ = 0, we 
obtain /' = i' s . 

The dynamic Neel texture leads to Gilbert damping in 
the CAF bulk. The amount of spin current lost in the 
bulk reads I s - /' = as£lLW flOlL where a is the bulk 
Gilbert damping parameter (whose microscopic origin is 
discussed at the end of the paper), and s = fiS/y is the 
saturated spin density, with Y denoting the area per spin 
of the CAF. The global frequency is then given by 

Ml = ---eV_ , (5) 

y + Y + y a 

where y a = 4nasLW/fi. Then the amount of spin cur¬ 
rent generated on the detection side by the superfluid spin 
transport reads 


where /m z = ~^m z F is the force conjugate to M z and F is 
the free energy of the CAF [in obtaining Eq. ([3]), we have 
assumed a symmetry 5? of the device in FigJTJa) under 
time-reversal followed by a n spatial rotation about the v 


I' = i- TL - eV _ . (6) 

4 ny + Y + y a 

Eqs. |5]l and ((6j> constitute the main results of this work. 
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FIG. 4. The line junction on the injection side. Charge currents 
entering vertices a and b redistribute according to scattering prob¬ 
ability matrix S . Mixing of charges between the two channels in 
the line junction is quantified by an effective conductance per unit 
length g(y). The green box denotes the injection region. 


This phenomenological result can be derived in a rotating 
frame, in which the spin spaces of all the edge electrons 
and the CAF rotate about the z axis with frequency £1 In 
this frame, the voltages of the edge channels emanating 
from the reservoirs are shifted as V a —> V a = V (T -crhQ./2e 
and V' (T —> V' T = V' r - crfiLl/2e. Since the Neel texture is 
static in the frame, the full spin current injected into the 
CAF on the injection side can be obtained by substituting 
Vo- in for V a in the expression for the static contribution, 
i.e., I s . With this substitution, one arrives directly at the 
result Eq. @041. 

Kinetic theory for injection/detection regions. —Here, 
we develop a simple microscopic model for the parame¬ 
ters y and y'. On the injection side, y quantifies the extent 
to which the two edge channels equilibrate inside the injec¬ 
tion region. Within linear-response, y can be evaluated for 
the (static) CAF in equilibrium. Due to the adjacent CAF 
order, the spin states along the line junction may deviate 
away from the ±z directions due to the effective field cre¬ 
ated by the CAF [see blow-up in Fig.[lja) and Fig.|3). The 
spin quantization axes there are thus expected to be canted 
away from the ±z axis and we label them ff and JJ,. At ver¬ 
tices a and b , the relative spin misalignment between the 
(t, i) and (ff, |l) states, together with sources of momen¬ 
tum non-conservation there, e.g., edge disorder and the 
sharp directional change of the edge, can give rise to inter¬ 
channel charge scattering. The redistribution of charges 
at these vertices must obey charge conservation, and can 
be parameterized by an energy-independent transmission 
probability t e [0,1] (under the assumed symmetry 2?, the 
two vertices are characterized by an identical probability) 


4 ( 0 ) \ 
4(0) ) 


- 8qS 




= S 


(hm\ 

\ h(W)) 


(7) 


where /^(y) (with <u =f[, JJ,) is the local charge current 
flowing along the line junction in edge channel <u, 5 = 
t &o + (1 - t)& x is the scattering probability matrix at the 
vertices, and o" 0 and fr x are the 2 x 2 identity matrix and 
the x component of the Pauli matrices, respectively. 

For inter-channel scattering inside the line junction, one 
requires: (i) spatial proximity of the two channels, such 
that there is sufficient overlap of their orbital wave func¬ 
tions; (ii) elastic impurities, providing the momentum non- 



p = W/L 

FIG. 5. Effective spin conductance G* ff = I'JeV- as a function 
of the aspect ratio p = W/L. We fix the edge equilibration length 
l and the system length L so that l = L/t - 10. The red and 
blue curves are, respectively, for full ( t - 1/2) and no (t = 1) 
inter-channel mixing at the vertices, and effective Gilbert damp¬ 
ing or e ff - 0-05 is used. The solid and dashed black lines, respec¬ 
tively, correspond to the weak (ar e ff = 0.05) and strong (or eff = 5) 
Gilbert damping with a transmission probability at the vertices of 
f = 0.9. 


conserving mechanism necessary to overcome the mis¬ 
match in Fermi momenta of the two channels; and (iii) 
a spin-flip mechanism, assumed here to be provided by 
the neighboring CAF. All three factors go into defining 
the tunneling conductance g(y) per unit length between the 
edge channels. In terms of g(y), the change in current on 
channel or is given by SI M (y) = ^giyW^iy) - V^(y)]6y, 
where V® is the local voltage on edge channel or [we as¬ 
sume that the edges are always locally equilibrated at all 
points y such that the voltage at each point is related to the 
local current through V^iy) = h(y)lgQ \• Then, the cur¬ 
rents inside the line junction satisfy 


dk 

dy 


d_h 

dy 


lk(y)-h(y)l 

Sq 


(8) 


Assuming a position-independent tunneling conductance g 
and defining the edge equilibration length £ = gg/2g, the 
currents entering vertex b is then given by 


I h(W) 
\ h(W) 


1 / 1 + e- w ' ( 1 - e~ w/e \ I 4(0) \ 

2 1 1 - e- W!t l+e- w ' e ]U(0) J ‘ 


(9) 


Combining Eqs. 0 and ([9]), the parameter y (on the injec¬ 
tion side) reads 


y = 1 - (1 - 2 tfe~ m . (10) 

A fully analogous consideration on the detection side leads 
to y = 1 —(1—2 t') 2 e~ w/e \ where t' is the transmission prob¬ 
ability at vertices a ' and b', and V is the edge equilibration 
length associated with the line junction on the detection 
side. 

Discussion .—For clarity, the results are now discussed 
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for the symmetric case, in which the injection and detec¬ 
tion sides are characterized by an identical transmission 
probability and edge equilibration length, i.e., t = t' and 
€ = €'. Recall that € describes the length scale over which 
the two relatively-biased edge channels in the line junc¬ 
tions chemically equilibrate via inter-channel scattering; 
we fix € and the length of the CAF region L to the ratio 
l = Ljt = 10. The effective spin conductance through the 
CAF, G* ff = I'JeV _ [see Eq. @], can then be expressed in 
terms of two variables: the aspect ratio p = W/L and the 
effective Gilbert damping ar e ff = asL 1 /Ti in the bulk CAF. 

The effective spin conductance is plotted as a function 
of the aspect ratio p for different t and a e ff. Full mixing of 
the edge channels at the vertices, i.e., t = 1/2, entails spin 
current injection only at vertex a. In this case, increasing 
p only increases the effects of Gilbert damping, since the 
latter is a bulk effect of the CAF, and the spin conductance 
monotonically decreases essentially as G* ff oc p~ l (see the 
red line). We call this “damping-dominated” behavior. If 
no scattering occurs at the vertices, i.e., t - 1 , spin current 
can only be injected within the line junction. For widths 
much smaller than the equilibration length, i.e., pi « 1, 
increasing the width gives an enhancement in the injected 
spin current that overcomes losses due to Gilbert damping, 
and a linear increase G s q& oc p (see the blue line) is ob¬ 
tained. However, as the width increases beyond the equi¬ 
libration length, spin injection no longer increases while 
Gilbert losses continue to increase. This leads to the even¬ 
tual decay G* ff oc p~ l for large p. We call this “weak- 
damping” behavior. 

For partial inter-channel mixing at the vertices, 0 < t < 
1, G^ ff has a qualitatively different dependence on p for 
different Gilbert damping strengths. Let us consider de¬ 
vice widths much less than the edge equilibration length, 
W € (the blue shaded region in Fig. [5]). Here, the 
Gilbert damping effects are small as long as the effective 
spin conductance in the lossless regime, characterized by 
y, is much larger than the effective spin conductance y a as¬ 
sociated with Gilbert damping. For y » y a , G^ ff exhibits 
the weak-damping behavior where an enhancement in the 
spin conductance is observed for W <$c t (see the solid 
black curve). Damping-dominated behavior is restored for 
stronger Gilbert damping, i.e., y <k y a , where the effective 
spin conductance monotonically decreases with the aspect 
ratio (see the dashed black curve). 

The Gilbert damping parameter a, as introduced in 
Eq. quantifies macroscopic relaxation of spin angular 
momentum polarized along the z axis, arising only in the 
presence of both spin-orbit interactions and microscopic 
degrees of freedom for energy dissipation. For a perfect 
graphene membrane, the spin-orbit coupling is known to 
be weak, while energy dissipation may be provided by 
phonons and magnons. Gilbert damping should vanish in 
this case as one approaches zero temperature, as phonons 
and magnons freeze out. Magnetic or heavy-element im¬ 


purities and/or enhanced spin-orbit interactions due to sub¬ 
strate can increase the spin relaxation rate. It is known, for 
instance, that spin-orbit interactions in graphene can be en¬ 
hanced, e.g., by hydrogenation m or interfacing it with a 
heavy-element-based semiconductor GE Additional dis¬ 
sipation channels, furthermore, can stem from substrate 
phonons or an ensemble of two-level systems (rooted in, 
e.g., crystalline defects either in graphene or the substrate). 
The latter could lead to energy relaxation (and thus finite 
a) even as T —> o m. In a flat graphene membrane, 
U(l) symmetry-breaking magnetocrystalline anisotropies, 
which, in principle, could pin the orientation of the in¬ 
plane staggered magnetization and thereby quench the su¬ 
perfluid state (9] 101, are expected to be small due to 
graphene’s weak spin-orbit coupling and high space-group 
symmetry of the hexagonal lattice. 
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